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ABSTRACT 

The eigenvalue problem for the linear stability of Couette 
flow between rotating concentric cylinders to axisymmetric 
disturbances is considered. It is shown by numerical 

caluculations and by formal perturbation methods that when the 
outer cylinder is at rest there exist complex eigenvalues 
corresponding to oscillatory damped disturbances. The structure 
of the first few eigenvalues in the spectrum is discussed. The 
results do not contradict the "principle of exchange of 
stabilities"; namely, for a fixed axial wavenumber the first mode 
to become unstable as the speed of the inner cylinder is 
increased is nonoscillatory as the stability boundary is crossed. 
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1. INTRODUCTION 


In this paper we reconsider the classical eigenvalue problem 
for the linear stability of Couette flow between rotating con- 
centric cylinder to axisymmetric disturbances. First we state the 
mathematical problem. 

Let r, 6, z denote the usual cylindrical coordinates, and 
le.t R_]_, ^2' ^2 radii and angular velocities of the 

inner and outer cylinders, respectively. In the formulation of 
the eigenvalue problem we will not take '^2 squal to zero; however, 
in all the calculations ^2 equal to zero. With d = R 2 - Ri» 
we define the following dimensionless parameters: 

n = Rx/R2 / ^ ~ d/R 2 , u = ^ 2 /^!' 

( 1 . 1 ) 

Y = 2n(l “ vi/n-)/(l + n ) , R = ^^iRl^/v, 

where v is the kinematic viscosity. We scale all length variables 
with respect to d, time with respect to d^/v, and all velocities 
with respect to Rx^^x* Then the linear stability problem for the 
stability of Couette flow to disturbances proportional to 
exp(at + iaz), wnere a is real and positive, can be written in 
the form 

(DD* - a^)2u - 2a2RFjj^(r)v - a(DD* - a^)u = 0 , (1.2) 

yRu + (DD* - a^)v - ov = 0 , (1.3) 

for rx < r < r 2 with 

u=Du=v=0atr= rx/ ^ 2 * (1«^) 
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Here 


D=d/dr, D* = D + r-1, ri = n/(l-n), r 2 = ri + 1; (1.5) 

u(r) and v(r) are proportional to the radial and azimuthal 
components of the disturbance velocity; and Fj^(r), the dimension- 
less Couette angular velocity, is 

• = ■ T + ~ ' (1.6) 

(l+n)(l-n)^ r^ 

Equations (1.2) - (1.4) define a non-self ad joint eigenvalue 
problem 

h(n,P,a,R,0)=O. (17) 

DiPrima and Habetler (1969) have shown that for fixed values of a, 
R, p, and n this eigenvalue problem has a countable spectrum {oj}, 
which can be ordered with Re( 0 i) > Re(a 2 ) > •••, with no cluster 

points in the complex plane and that the corresponding 
(generalized) eigenfunctions span a certain Hilbert space. 

In this paper we will show by numerical calculations and by 
formal peturbation methods that for fixed p and p = 0 (outer 
cylinder at rest) there are values of a > 0 and R > 0 such tnat 
some of the 0 j are complex. To the best of our knowledge, this is 
the first demonstration of the existence of complex eigenvalues 
for the boundary value problem defined by Eqs. (1.2) - (1.4). This 
result contradicts a proof by Yih (Main Theorem p. 299, 1972b) that 
all the eigenvalues are real when the cylinders rotate in the same 
direction and the circulation of the basic flow decreases in the 
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outward radial direction. However, all of the complex eigenvalues 
we have found correspond to damped disturbances, Re (a) < 0. Thus 
the results do not contradict the conjectured (and widely 
believed) principle of exchange of stabilities; namely, for fixed 
values of n, a > 0, and u > 0 the first eigenvalue (oi) to cross 
the imaginary axis as R is increased is real. However, the present 
results show that a proof of- the principle of exchange of 
sta'bilities must be restricted to a study of the behaviour of the 
first eigenvalue. The condition that cr^ is real and simple is 
required in a rigorous proof of the existence of Taylor vortex 
flow following the instability of Couette flow; see Velte (1966) 
and Kirchgassner and Sorger (1969). 

In Section 2 we give an example of how complex eigenvalues 
can arise due to small perturbations of a selfadjoint eigenvalue 
problem. The example is related to the eigenvalue problem (1.2)- 
(1.4), but is simple enough that one can carry out the 
calculations readily and explicitly. In Section 3 we give some 
results for the numerical calculation of complex eigenvalues of 
Eqs. (1.2)-(1.4) and confirm these results by formal perturoation 
calculations. In Section 4 we study tne eigenvalue problem in the 
limit a > 0 in order to obtain a better understanding of how the 
complex eigenvalues occur. Finally, the results are discussed in 
Section 5 where we also make a few remarks aoout the limiting case 

a 

For many purposes, and especially for numerical work, it is 
convenient to write Eqs. (1.2) and (1.3) as a system of six first 
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order equations. If we introduce the axial component of the 
perturbation w = ia“-^D u and the pressure perturbation p and let 
q = [p, dv/dr, dw/dr, u, v, w], we can write Eqs. (1.2) and 
(1.3) in the form 


dq/dr - Aq - oBq = 0, 


ri < r < r 2 


( 1 . 8 ) 


with the boundary conditions 


u = v= w = 0 atr = ri, r 2 < 


(1.9) 


The matrices A and B are 


and 


A= 


B = 


0 

0 

ia 

0 

0 

.0 

0 

0 

0 

0 

0 

0 


0 

-i/r 

0 

0 

1 

0 

0 

0 

0 

0 

0 

0 


-la 

0 

- 1 /r 

0 

0 

1 

0 

0 

0 

0 

0 

0 


-.2 


a' 
Ry 
0 

1 /r 

0 

0 


-1 

0 

0 

0 

0 

0 


0 

1 

0 

0 

0 

0 


2RF^(r) 

a^+l/r^ 

0 

0 

0 

0 


O' 

0 

1 

0 

0 

0 ; 


“\ 

0 

,2 


-xa 

0 

o/ 


( 1 . 10 ) 


Finally, for use later, it is helpful to have the results 


Y = (1 - u) - 6 ^ - 52 1 + 7u 4 . 0(^2) 

2 2 


( 1 . 11 ) 
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F,(r) = 


and, for u = 
2YF^(r) 

where x = r 


4 [a . - (1 - u)x] . . 3 U - u)x 

2 2^2 

2 3 

(1 - y)x ] + 0(6 ) ' ' 

0 , 

= 6(1 - X) - 6^(x - x^) + 0(6^) , ( 

-rj^ = r + 1-1/6. 


. 12 ) 


.13) 
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2. AN EXAMPLE 

If we solve Eq. (i.3) for u, substitute in Eq. (1.2), let x = 
^ ~ ^1' “ R^6(l - u^) f and take the small-gap limit 6 0 with 

o, T, and a fixed we obtain the classic small— gap equation 

(d2/dx2 - a2 - a)2 (d2/dx2 - a,2 

+ a2T (1 - £x)v =0, 0 < X < 1, 

where e = (1 - + y ) . We consider Eq. (2.1) with the 

boundary conditions 

V = v'' = v'*'' = 0 at X = 0 and x = 1. {2.2) 

The eigenvalue problem defined by Eqs . (2.1) and (2.2) is the 
small-gap limit of the auxiliary problem introduced by Yin [eqs. 
(11), (12), and (14), 1972a] for a > 0 , t > 0 , 0 < e < 1 . Yih (o. 

296, 1972b) asserts that for > 0 (T > 0), a > U, and 0 < e < l 
the eigenvalues a of his auxiliary system are real. We will show 
that this is not true for the corresponding small-gap equations 
(2.1) and (2.2). 

For e = 0, Eqs. (2.1) and (2.2) are simply those of the 
classic Benard problem with free-free boundaries - a selfadjoinr 
boundary value problem. The eigenvalues a can be split into two 
sets and where 


= 


2 2 
_ n 7T 


- a 


* n^Ti^ + 


n = 1, 2, 


(2.3) 
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The corresponding eigenfunction for and is Vj^(x) = 

sin mix. Notice that if T = 0, then the eigenvalues are negative 
and each eigenvalue has multiplicity two ( 0 ^'*’ = On~) with only one 
eigenfunction; however, this degeneracy is not of interest at tne 
moment. More important is to observe that for any n, 0 ^'*’ 
increases monotonically with increasing T and eventually becomes 
positive while decreases montonically with increasing T. 

Thus we can choose positive integers N and M with N > M such 

that 




= - a2 




2 

a Tn 

N2n2+a^ 


(2.4) 


where 


2 2 2 
^ (N^ - 

a 


1 + 1 

(n2tt 2 + a2)i/2 (m2tt2 + a2)V2 


-1 


.(2.5) 


For a given value of a with e = 0 and T = Tq the eigenvalue 
problem (2.1) and (2.2) has an eigenvalue of multiplicity two, a = 
with two linearly independent eigenvectors vi>j(x) = 
sin Nttx and vm(x) = sin Mux, N t M. We will now show by standard 
perturbation methods that for values of T close to Tq and e 0, 
Eqs. (2.1) and (2.2) have complex eigenvalues. 

We write 


a = a, 


+ + 


= -(N^tt^ + a^) + 




2 

a Tn 

r2,2 


T 


+ eTj^ + • • * , 


( 2 . 6 ) 


v(x) = Vq(x) + evj^(x) + ••• 


/ 
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where Vq(x) - a sin Njtx + 0 sin Mux = av^(x) + 0v.j(x). The 

constants a and 0 are to be determined as part of the calculation. 
Substituting the series (2.6) in Eqs. (2.1) and (2.2), we obtain 

Lvq = [(dVdx^ - a^ - Oq)^ (dVdx^ - a^) 

^ '"O " ° ' (2.7) 

and 

LVf = 2oj_(d2/dx2 - a2)2v^ - 2oQaj^ (d^/dx^ - h2)Vq - 

+ a^T^xv^ , (2.8) 

and the boundary conditions (2.2) for Vq and Vj^. Equation (2.7) 

with the boundary conditions (2.2) is automatically satisfied. In 
order for tne boundary value problem for v^^ to have a solution, it 

is necessary that the right side of Eq. (2.8) be orthogonal to 
vn(x) and Vv^(x); see Courant and Hilbert (pp. 346-350, 1953) or 
Case I of the Appendix of this paper. This leads to two linear 
homogenous equations for a and 0, and the condition that these 
equations have a non-trivial solution is 
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It is clear from Eq. (2.9) that if N + M is odd and 
(Ti - Tq/ 2)2 < 64N2M2To2/ir - M^) then is pure imaginary. 
Thus for some, but not all, values of Tj_ the perturbation in e of 
the double eigenvalue, gives complex eigenvalues. For example, if 
Ti = Tq/ 2, then 



4NMaTo ^ 

tt2(n^ - m2)^(N^tt2 + a2)-L/4 (m2,j2 + ^ 2 ^ 1/4 


( 2 . 10 ) 


Hence for e 0 with T = Tod + e/2) there are complex eigenvalues 
of the boundary value problem (2.1) and (2.2) of the form ± i£v 

+ 0(e2). While the above arguments do not provide a rigorous 
proof that the boundary value problem (2.1) and (2.2) has complex 
eigenvalues for certain positive values of a, T, and e it is 
strongly suggestive of such an assertion. Moreover, it suggests 
that the corresponding problem without the small-gap limit, tne 
auxiliary problem discussed by Yih, also has complex eigenvalues 
which contradicts his assertion. 
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3. COMPLEX EIGENVALUES 

Vve have solved the eigenvalue problem (1.8), (1.9) for 
several different sets of parameter values using a standard 
shooting procedure. The integration routine was a fourth order 
Runge-Kutta procedure with a step size of 0.025, and Mueller's 
method was used for the eigenvalue search routine. In Figure 1 
we. show our calculation of the first five eigenvalues for n = 
0.877, li = 0, and R = 150. For this value of n the critical value 
of R, above which there exist values of a such that a > 0, is R^ = 
119.3. The corresponding critical value of a is a = 3.13. 

It can be seen in Figure 1 that the first eigenvalue is real 
and becomes positive for a finite band of wavenumbers 
(approximately 1.7 < a < 5.3) corresponding to the band of 
wavenumbers inside the neutral curve at R = 150. The second and 
third eigenvalues are real and negative for R = 150 and 1 < a 
< 8. The fourth and fifth eigenvalues merge at a = 2.56 and form 
a pair of complex conjugate eigenvalues for 2.56 < a < 6.6. 

Indeed for a < 1 the second and third eigenvalues cross over and 
for a very small interval of wavenumbers form a pair of complex 
conjugate eigenvalues. In Figures 2 and 3 we show the first five 
eigenvalues for n = 0.75 and n = 0.5, respectively, for R = 150. 
The overall structure of the eigenvalues remains virtually 
unchanged as n is varied, except that the interval of wavenumbers 
for which the fourth and fifth eigenvalues for.m a complex 
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conjugate pair decreases with decreasing values of n. We found 
no indication of a possible crossing of any eigenvalue with the 
first eigenvalue in any of our calculations; note that R = 150 
corresponds to 1.26R(^r 1.75 Rq, and 2.20R^ for n equal to 0.877, 
0.75, and 0.5, respectively. 

At a value of a, say ao^ in Figures 1, 2, and 3 where two 
real eigenvalues merge to give an eigenvalue of multiplicity two 
our calculations show that there is one only eigenvector. A 
generalized eigenvection can be obtained by numerical 
integration of the equation 

d(j;/dr - Aij; - aBi|; = Bq (3.1) 

with the boundary conditions (1.9) and where q is the eigenvector. 
As a check on the consistency of the numerical calculations we can 
follow tne procedure described in the Appendix (Case II) to cal- 
culate the variation of the two eigenvalues with a for [a - aol 
small. In addition to calculating q and we must also cal- 
culate the eigenvector and the generalized eigenvector of the 
adjoint system 

dq*/dr + Atrq* + =0 (3.2) 

with the boundary conditions that the first, second, and third 
components of q* vanish at r = rq and r = r 2 » 

We find for n = 0.377, y = 0, R = 150, and ao = 2.5576 thac 


04^5 = -92.1318 ± 3.229 (aQ - a)V2 + o(a - ag) . 


(3.3) 
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Formula (3.3) shows that for jag - a| small, the eigenvalues 
and are real for a < a^ and are complex conjugates for a > a^. 

The agreement of the perturbation formula (3.3) with the numerical 
calculations is satisfactory, especially for a^f as can be seen in 
Figure 4. 

It is impossible to carry out a complete analysis of the 
spectrum of the eigenvalue problem defined by Eqs. (1.2) - (1.4) 
or Eqs. (1.8) - (1.9) for all allowable values of n, u, a, and R. 
Our primary goal was to show the existence of complex eigenvalues 
for VI = 0 with values of a and R typical of those for which there 
exist eigenvalues which are real and positive (unstable modes). 
Having found complex eigenvalues we snail now consider the 
limiting case a ->■ 0- We will find that the structure of the 
spectrum found by numerical calculations for the full equations is 
preserved in this limit. As a consequence, we can gain some 
understanding of the origin of the complex eigenvalues. 
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4. THE STRUCTURE OF THE SPECTRUM IN THE LIMIT a > 0 

In this section we consider the eigenvalue problem (1.2) 

- (1.4) in the limit a -► 0. It is Icnown that on the neutral curve 
R = 0(l/a) as a ->■ 0. If we take R = 0(l/a) then it is also clear 
from Eqs. (1.2) and (1.3) that in order to obtain a meaningful 
problem we must take u = 0(a). It is also convenient to introduce 
the scaling that is used for the small-gap problem. Thus we let 

R^ = aR, u = aU, V = yRa^. (4.1) 

If we substitute these expressions in Eqs. (1.2) - (1.4), and then 
let a -»■ 0 with n , o, Ra/ H/ and V fixed, we obtain the eigenvalue 
problem 



with 

U=DU=V=0, (4.3) 

and where G£(r) = 2Y6“^F£(r). 

We will first study the eigenvalue problem (4.2) and (4.3) 
for small values of 6 and Ra^iS using perturbation methods. For 
this purpose we set x = r - r^/ observe from Eq. (1.13) that G^^(r) 
= (1 - x) + 6(x - x^) + O(o^), and then let o 0 and Ra^5 -<• 0 in 
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Eq. (4.2). We obtain 


'd^/d; 


d2/dx^ 



dVdx2 0\ / u 


with the boundary conditions 

U = dU/dx =V=0 atx=0 and x = 1. (4 

This is a non-self ad joint eigenvalue problem with the following 
sets of real eigenvalues and eigenvectors: 


'"An = = 1 - cos 2n:rx, = 1_^ (1-cos 2n7rx) 


4 n^ TT^ 


- - sin IniTx ' 

4n7T 


“Bn “ ■ " ' ' “Bn “ ''sn “ ' 


“Cn “cn ' ^ - 


_ Sin An / , 

^ H_ (A x - sin Ax, 

i-COS An ^ 


n = 1, 2, .... The A^ are the roots of A/2 = tan A/2 witn A^^ = 

8.936, A^ = 15.45, and A^ (2n + I):: for n ->• ®. The functions 

V are the solutions of 
Cn 


(dVdX^ ^ ^ 


(4.5) 
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The first nine eigenvalues are 


<^1 

^5 

^6 


®B1 

= -IT = - 

-9.870 

®A1 

(II 

CN 

1 

II 

-39.48 

*^B2 

= -4.2 H 

-39.48 

'^Cl 

= -X^2 , 

-80.75 


Og3 = -9 ti^ = -88.83 
a,T = -16n^ = -157.9 


= -16n^ = -157.9 

^8 ^ ^ C 2 " ~^2 ~ 

ag = Ogg = -25 tt^ = -246.7 


The reason we have labeled the eigenvalues in the above manner 
will become clear later in this section. For the moment, notice 
that with the exception of the first eigenvalue, the B eigenvalues 
for n even pair with the A eigenvalues and for n odd pair with 
the C eigenvalues. 

It is possible to study the behavior of the multiple 
eigenvalues at a = — 4m^ir^, m = 1, 2, ••*, for small values of Ra^ 6 
and 6 using formal perturbation methods. In order to study the 
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effect on these eigenvalues of comparable perturbations in Ra 2 5 
and (5, we set 

= a2s2, A = 0(1) (4.10) 

t.n Eq. (4.2), and then expand in 6 following the procedure 
described in Case I of the Appendix. In order to carry out this 
calculation we need the adjoint eigenvectors for the eigenvalue 
problem ( 4 . 2 ) and ( 4. 3 ) . If we let q = [u, V]tr ^„d introduce the 
inner product 


< q 


1 ' 


> = 


(^ 1^2 + Vj^V^) dx. 


(4.11) 


tnen we find that the adjoint eigenvalue problem is given bv Eq. 
(4.4), with each matrix replaced by its transpose and Eq. (4.5). 

The adjoint eigenvectors corresponding to a, = fu , V and 

tr 


■1 


Bn 



V 


Sn- 


are 


An 


= 1 - cos 2mrx, V 


An 


= 0 , 


n — 1 , 2 , 


(4.12) 


and 


U 


Bn 


= 1 - cos 


4n 2 ti 3 


hit X _ sin mix ^ x cos mrx 


2n‘^TT^ 


2n 


T,3 


Bn 


= sin nTx, 


n = 1, 3, 


(4.13) 


U 


Bn 


★ 


3 

X 

2n^ IT 3 


(cos 2n7Tx - 1) , 



sin nnx. 


n = 2, 4, •• 
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If we let 


P = 


dVd: 


0 


1 


then one can readily verify the orthogonality relations 


(4.14) 


and (4.15) 

Returning to the eigenvalue problem (4.2) and (4.3) we now 
consider the perturbation of the double eigenvalue a = -4m-Tt for 
small 0 witn A =0(1). The series expansions are 


= 4m + 5 u + 6 + 


(4.16) 


. 2 , 


q = + s qe,2m> + «qi + « qj + ' 


for m = 1,.. 2, •••, where the ratio a/3, which depends on m, is to 
be determined in the course of the analysis. (vie use u and v as 
perturbation coefficients. There should be no confusion with 
their earlier usage for the ratio n2/^i' kinematic 

viscosity. ) Following the procedure described in Case I of the 
Appendix, we find at 0(5) that u = 0 and 


“ 


X (cos 2 raTTx - 1) 

2 


Sm^ir 2 


- 2Lf_ sin 2mTrx 
8m u 


(1 - cos 2mTTx) 
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+ B 


_ 1 


± X sxn 
2 


2imr X 


(4.17) 


The function is only determined up to additive multiples of the 
eigenfunctions and explained in the Appendix 

there multiples can be chosen to be zero. If they are included in 

qj^,. they will not affect the calculations of v. At 0( 6^) we 
obtain the following two linear homogeneous equations for a and S: 


a 


_J + A_^ / 1 + 105 

32mTr 256m5T;^ ISm^TT^^ 



3m^Tr^ _ 3A^ _ 2m^T\^v 

^ 16m^TT 2 


+ B 


- 3jLl 

8mTT 


= 0- 


0» 


(4.13) 


These equations determine v and the ratio a/B. There are two 
values of v given by 


/-3a2 * 

’ 9 _ 

27a2 

+ 6a4 I 

5 + 105 1 1 

V2. 

\ 32m'iTr4 

¥ 

64m ^ 

1024m^Tf8 [ 

^ 16m2^2|j 

i 


(4.19) 


For 5 0 with 5 >> whicn corresponds to letting A 0, 

we find that v = ± 3/4. Thus a perturbation in o splits the 
double eigenvalue at a = -4m^u^ into two real and unequal 
eigenvalues , 
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= — ^ + 0(A^)J. 


(4.20) 


For Ra + 0 with Ra^ >> 6, which corresponds to letting A -»■ ®, 
we find that 


= A- 


6 4m^TT 


- 3 ±Vs I + 1£5 


^ IGm^ir ^ j 


1/2 


+ 0(1/A^) 


Hence, again the double eigenvalue at a = -4m^u2 is split into two 
real eigenvalues, 


.2.2 , 2_ 1 


a = -4n; TT + (R 6 ) 

^ 64m 2 




5 + 105 


1/2 


2 16m^ ir2 


+ 0(1/y ) 


(4.2i: 


However, for each value of m there is a finite interval of values 
of a 2 given by 


5 + 105 


m^TT 2 


-72 


ISm^TT 


/ 




+ 384 < 0 


(4.22) 


for which Eq. (4.19) yields complex conjugate values for v. For 
m = 1, the values of v are complex conjugates if 


8.54tt4 < a 2 < 14.21:i-^. 


(4.23: 


We can interpret this result as follows. Suppose we fix 5, 
0 < 6 <<1, and take Ra = 0 so Ra^o = 0 and hence A = 0. 
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Then the double eigenvalue a = -4m^Ti2 of Eqs. (4.2) and (4.3) when 
6=0 splits into two real eigenvalues, a = ± 352/4 at 

leading order. Now we increase but with << 62 and hence 

A << 1; then the dominant correction term to a = -4m2ir2 remains 
± 36 2/4 and the two eigenvalues are still real. However, as Ra is 
increased still further until Ra^6 = 0(62) the correction term v 
becomes complex when A reaches the lower value determined by Eq. 

(4*. 23) and v stays complex until A reaches the larger value 
determined by Eq. (4.23). For 62 << « 1 so a 2 >> 1 the 

eigenvalues are real and distinct differing from -4m2ir2 by terms 
0(Ra^6) as given by Eq. (4.21). 

In Figure 5 we show the results of the perturbation analysis 

for 6 = 0.123 (n = 0-877) for the eigenvalues 02 and 

corresponding to m = 1. For R =0 the eigenvalues are a-, = 

-39.4897 (corresponding to ) and = -39.4671 (corresponding 

2 

to Om ) ; as R is increased tney coalesce at R = 10.10 (A = 
0^3 3 

8.54tt4) and become complex conjugates until R^ = 13.05 ( a2 = 
14.21ir4) ait which point they coalesce again; and for R^ > 13.05 
they are again real. Also in Figure 5 we show the corresponding 
results obtained by solving the full equations (4.2) and (4.3) 
numerically. We have labeled the modes according to their 
ordering at 6 = 0 and R^ = 0. Note the crossing of the 
eigenvalues. The perturbation and numerical results are in 
reasonable agreement even though R^^o is not numerically small for 
a2 = 8.54tt2 and larger values of a2. We note that it would not 
have been easy to determine these eigenvalue curves numerically 
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without the information gained from the perturbation analysis. 

For the other doable eigenvalues a = -4m^Tr2, m = 2, 3, at 

6=0 and Ra ~ ^ splitting structure for perturbations in 6 

and Ra is similar to the case m = 1 just discussed. 

For a fixed value of 6 and moderate values of Ra it is 
necessary to resort to standard numerical procedures (discussed in 
Section 3) to solve the eigenvalue problem (4.2) and (4.3) for the 
variation of the eigenvalues a with Ra« Iri Figure 6 we show the 
development of the first five eigenvalues with Ra for 0 < Rg 
< 1200 for n = 0.877 (6 = 0.123). The modes are labeled 

accordingly to their ordering at. Ra = 0* Along the branch denoted 

* - ** the eigenvalues and are complex conjugates. The 
corresponding modes interchange their order at the point ** just 

as was the case for the 02 eigenvalues. Our numerical 

calculations (and perturbation calculations for 02 and ) show 
the following continuous mode association: ^ 2 '^ *^A1 ' 

° 3 % 2 ' ^ 4 ^ *^ 01 ' % 3 ‘ 

For the modes we have investigated, we find tnat the 3 
eigenvalues initially increase monotonically witn Ra/ while the A 
and C eigenvalues initially decrease with Ra* When a B and an A 
eigenvalue or a H and a C eigenvalue coalesce (recall the pairing 
mentioned earlier) the eigenvalues become complex conjugates for 
some finite interval of values of Ra/ the eigenvalues then 
coalesce again and split as two real eigenvalues with the 
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corresponding modes interchanging their order. The only 
eigenvalue which does not have such a crossing (and hence for some 

Ra an association with a complex eigenvalue) is txie first, or most 

unstable mode, None of our numerical calculations 

showed the existence of a complex eigenvalue a = o + ia , 

r i' 

with Oj- > 0; however we did not make an exhaustive search. 

Finally, we want to mention a not so obvious way in which the 
B and A modes can cross. This is associated with a branching at 
points at which da/dR^ = 0. At such points the inverse eigenvalue 
relation R = R(a) has a double eigenvalue which can split apart as 
a pair of complex conjugate eigenvalues as c (real) is varied. 

Thus in Figure 6 the points Q and O at which da/dR = 0 are 
connected by a branch on which a is real and there are two complex 
conjugate values of R^. In this way the A1 and the 33 modes 
interchange their order with a remaining real but R^ being 
complex. Of course for the physical problem R, and hence R^, is 
real . 

Corresponding to the segment * - ** in Figures 5 and 6 the 

eigenvalues 02 and and and a^, respectively, are complex 
conjugates. From the numerical calculations these intervals 
correspond approximately to 


and / 

« 

CO 

< R, < 9.8 
a 

or 

7 . 8 /a 

< R < 9.8/a; 

and , 

320 

< R < 720 
a 

or 

320/a 

< R < 720/a. 


With this in mind we investigated the full problem for n = 0.877 
to determine regions in the a-R plane where complex eigenvalues 
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would be found. In Figure 7 we show this region (between the two 

curves) for the eigenvalues and The asymptotic results for 

a > 0 are in reasonable agreement with the numerical calculations 

even for moderate values of a. There is a similar region for the 

eigenvalues 02 and a^r and, we believe, for the other B 
eigenvalues as they cross the A or C eigenvalues . Again we note 

that the numerical calculations were motivated by the results of 
the asymptotic analysis. 
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5 . DISCUSSION 

We have demonstrated by numerical computation and formal 
perturbation methods that the eigenvalue problem for the linear 
stability of Couette flow to axisymmetric disturbances (the Tavlor 
problem) has complex eigenvalues when the outer cylinder is at 
rest. However, it is not inherent in the analysis that the outer 
cylinder is at rest (for example the asymptotic analysis for a 0 
in Section 4 ) and we anticipate the occurence of complex 
eigenvalues when the outer cylinder also rotates. All of the 
complex eigenvalues which we have found correspond to damped 
oscillatory modes; and in particular there is no indication that 
the first eigenvalue is complex for any values of n, a, and R with 
the outer cylinder at rest. Thus it seems likely that the 
principle of exchange of stabilities holds. However, the 
existence of complex eigenvalues shows that, unlike the situation 
for the Benard problem, a proof of this principle cannot be 
constructed by showing that all the eigenvalues are real. 

It is' interesting to speculate on whether an experimental 
investigation of the Taylor problem could confirm the existence of 
decaying oscillatory axisymmetric disturbances. In any 
experimental situation there are small imperfections and 
random fluctuations present which, presumably, will excite all the 
modes of linear thoery. However, all of the oscillatory decaying 
modes have large decay rates. The least damped oscillatory modes 

have h -4;t (corresponding to and 03 coalescing), but this 
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occurs for only a small area in the a-R plane. The eigenvalues 

and can be complex conjugates for a significant set of 
values of a and R, as is shown in Figure 7. For n = 0.877, a = 5, 

and R = 150 we find the complex conjugate eigenvalues* a = -110.2 
±i3.72. The physical time T for such a mode to decay to e~^ 
times its initial amplitude is given by T = 0.19/Oj^, which may 

allow sufficient time for observation. 

We have seen that the structure of the spectrum of the linear 
problem for a = 0(1) is preserved in the limiting problem 
discussed in Section 4 for a -<■ 0. Moreover, for the latter 
problem we found that the complex eigenvalues occur because there 


are two 

denumerable 

sets 

of 

negative eigenvalues at 

Ra 

= 0 which 

separate 

each other 

and 

move 

in 

opposite directions 

as 

R^, and 

hence R, 

is increased. 

One 

set 

of eigenvalues have 

Or 

increasing , 


and ultimately correspond to unstable modes, while Vj- decreases 
for the other set. It follows that there will be an infinite 
number of intersections and at these points of intersection we can 
expect the, occurence of complex conjugate eigenvalues. As we have 
noted these points of intersection correspond to damped 
oscillatory disturbances. Ic is of course possible that the 
eigenvalues corresponding to two unstable modes can grow at 
different rates with R, and hence eventually intersect so as to 
give rise to complex conjugate eigenvalues and growing oscillatory 


* This calculation was confirmed by Dr. P.M. Eagles using an 


independent program. 
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modes; however, we have not explored such possibilities. 

It is also possible to study the full eigenvalue problem, 

Eqs. (1.2) and (1.3), in the limit a « with n and R fixed. In 
this limit we can expect that a = 0(-a2) for all the modes so an 
expansion will take the form 

0 — a2(— 1 + a~^ y + ...)• 

It turns out that the eigenfunctions have a boundary layer 
behaviour at the cylinder walls and satisfy a fourth order 
differential equation in the interior. The parameter y and the 
higher order coefficients are determined by solving a sequence of 
fourth oraer equations with boundary conditions at the cylinder 
walls determined by matching conditions. We did not study this 
problem in detail, but our few calculations for the small-gap case 
did show that y can be complex and that the curves of Figure 7 
will asymptote to horizontal straight lines as a -► ®. It is 
interesting to note that if a ->■ <= with R = 0(a2), corresponding to 
the behaviour of R on the neutral curve, the asymptotic analysis 
of Hall (1982) shows that all the eigenvalues a are real. This 
result gives additional credence to the principle of exchange of 
stabilities, but a proof is yet to be constructed. 
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APPENDIX 

Consider a boundary value problem consisting of a system of 
ordinary differential equation of the form 

(L + eM + e 2 N + ♦••) q - o’ (A + eB + + •••) q = 0 , (Al) 

where L, M, N, A, B, C, are matrix differential 

operators with non-singular coefficients; q is a vector function 
satisfying separated homogeneous boundary conditions at the end 
points of some finite interval; say 0 < x < 1 ; e is a small 
non-negative parameter; a is the eigenvalue; and the order of the 
system (Al) is equal to the order of the operator L. The boundary 
value problem (Al), including the boundary conditions, need not be 
self ad joint. Our discussion is formal; moreover, we do not 
consider all possible cases. 

We assume that the boundary value problem 

(L - ah) q = 0 (A 2 ) 

has an infinite denumerable set of eigenvalue with oj_ = 02 

while all the other eigenvalues are simple. We are interested in 

how the corresponding eigenvalues 0^^ and 02 of (Al), whose limits 
are 0j^ and 02 as e -*■ 0 , vary with e. We consider two cases 
I. Corresponding to = 02 there are two linearly 
independent vectors q^ and q^ of (A 2 ), 

and 

II. Corresponding to Oj^ = 02 there is one eigenvector q^^ of 
(A 2 ) and a generalized eigenvector q^^j^ satisfying 
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(L - Oj_A) = Aqj^* (A3) 

We assume that qj^, q 2 / q 3 ' ••• in Case I and qj^, ^ 11 ' ^ 3 ' ••• 

in Case II are complete in an appropriate function space. Let 
<• , •> denote the inner product, and let 

(L* - ah*) q* = 0 (A4) 

witti appropriate boundary conditions be the adjoint boundary value 

problem corresponding to (A2). We also assume that ^ ^ 

small and non-zero. 

Case I . For this case the discussion is similar to that 
given in Courant and Hilbert (1953, pp. 346-350) for a multiple 
eigenvalue of a second order selfadjoint boundary value problem. 

We assume that the eigenvectors and adjoint eigenvectors have been 
chosen so that 

iff "ff 'if 

<Aq., <3j > = ^ ^ " hj * 

As part of the perturbation calculation it is necessary to 

determine the "directions" of the eigenvectors qj^ and q 2 of (Al) 
as e * 0. Thus we choose as our eigenvectors of (A2) the set 

^1 ^ “ll^l “ 12^2 ' *^2 "" “ 21^1 “ 22^2 /• 

Uj = q j , j = 3 , 4 , . . . , j 

where the a's are to be chosen so that q^^ + Uj^, q 2 U 2 as £ ■*• 0 . 
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We expand a^i <3n series 

On = cjn + ePn + + '•*/ + e^Wj, + 

n = 1 , 2, • ‘ . 

If we substitute (A7) in (Al), we obtain 
(L - anA)Vn = (“fi + On® ■*■ 

(L - OnA)Wn = (”N + a^C + w^® 

+ ( — M + a + y A ) V • 

Next we assume that we can write 


(A7) 


(A8) 


(A9) 


vn = 




Wj^ - 


j=l 


j=i 


^nj qj 


(AiU) 


If we substitute for (A3) and take tne inner product wicn 
q_j*, we obtain 

anm ( Ojn “ o^) = <(-M + ®n®^^n' ■*■ ^^n '^nm* 

For n = 1, 2 and m = 1, 2 the left side of (All) vanishes since 
0 l = 0 ^ and we obtain the equation 

i * * 

<(-M + (jj_B)qj_, qj_ > + u <(-M + a^B)q2/ ^ 

<(-M + aj^B)qj^, q2 > <("‘'^ + o q^ > + U 



-30- 


for U 2 « f'or n-m-3, 4, •••/ (All) gives u = 

^ n 

<(M - anB)qn' only need the first order 

corrections to the eigenvalues 7^ we can stop at this point. 
However, if we also require the first order corrections to the 
eigenvectors and the second order corrections to the 
eigenvalues, it is advantageous to reorganize our computations. 
We also find from (All) that 


- °|B)q,, q, > - 1 , 

- »iB)q2 ' 'Jl’> 

- g,B)q^ , q^*> 


= - cT^B)q^, q^% 

“22 <(M - aj^B)qj^, q^^*> - 

^ _ <(M - a^B)q q > _ ^ 


<(M - a^B)q^, q^ > 

where we have assumed <(H - Oj_B)q 2 , q^_*> ^ 0 and <(M - c B)q , 

* 12 

^2 ^ 0 • Now introduce new ad^joint eigenvectors u — 'I * ~ 

Sq 2 , u^ = -aq^ ^ ^2 ' ^ f or j = 3 , 4 , • • • , and choose 


*11 “22 “ 


As a result we he 


<Au^, u^ > = <u . , A u . > = 6^ . . 


We can also writs 


” oo 

Yj ^ Y ’ 


(A16 ) 
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Substituting this expansion for in (A8) and taking the 
inner product with Um*, we obtain 

” “ ^nm ^‘n "^nm (Al7 ) 

where = <(M - Oj^B)un/ will restrict our attention to 

the eigenvalues and 02 and the corresponding eigenfunctions; 
that is, n = 1 and n = 2. For n = 3, 4, ... the calculation is 

straightforward. Equation (A17) yields the following information 
for n = 1, 2 and m = 1, 2 


= F 


11 ‘ 


= F 


22 


12 


= 0 , 


21 


= 0 , 


(A18) 


and 


F ^ 1/2/ ® 3/4, •••. (A19) 


Equations (A18) are consistent with our choice of Uj^, 1 ^ 2 ' *^1' ^ 2 ' 
* ★ 

Uj^ , and U 2 as can be verified by doing the necessary algebra. 
Alternatively, if we had chosen the Uj according to (A6), the 

Uj* so that the biorthogonality condition (A15) is satisfied, and 
then used the expansions (A16 ) we would have obtained (A13) for 


and the ratios “21^‘^22‘ 

The coefficients Cj^^^ and C 2 j^ in tne series expansions for 
and V 2 are given by (A19) for m = 3, 4, We must determine 


*^ 11 ' '^ 12 ' ^^ 21 ' ^^22 


Imposition of the normalization condition 


<Aqn/ u^ > - 1 


n = 1, 2, 


(A20) 


yields c^^ = C 22 
(A9) for Wj^. 


0. To determine Cj ^2 *^21 


we must consider 
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Substitating the expansion for given in (Al6) and (a9) and 
taking the inner product witn Ujjj* gives 


CD 

j=l 


(A21) 


For n - m = 1 and n = m = 2 we obtain 


vn = <(N - 


anC 


nB.Un- C>- 

j=3 


- u 


(A22) 


and for n = i, ra = 2 and n = 2, m = 1 we obtain 


c 

nra 



<(N - o^C - y^BjUnf Ufn*> 


- V' ^nj ^ 3in 
3=3 


IA23) 


This completes the calculation of and 02 througn terms 

2 — — 

0( £ ) and the corresponding eigenfunctions and through terms 

0(e). To compute corrections to the other eigenvalues and 
eigenvectors as well as higher order corrections for a,, a,, q , 
and q^ is straignt forward in principle. In practice the 
calculations would probably not be carried out using eigenfunction 

expansions, but rather the nonhomogeneous equations for v^^ and v^ 
would be solved numerically after the solvability conditions (A18 ) 


had been used. 
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Case II For this case the discussions by Wilkinson (1965, 

pp. 62-70) for the matrix eigenvalue problem, Friedman (1956, pp. 

110-113 and 131-133), and Kato (1980, Chapter 2) are useful. Let 
★ ★ * 

qi , , ... and q^^j^ be the corresponding eigenvectors and 

generalized eigenvector of (A4). We assume that the vectors have 

been normalized (see Friedman for the case A is the identity 
operator) so that 


<Aqj_, 

* 

★ 

= <Aqii/ q^i > = 0, 

* ★ 
q^l > = <AqjLj^, q^^ > = 1, 


* 

•k 


<Aqj_, 

> 

= <Aqii» qj > = 0/ 

LJ. 

II 

U> 

• 

• 

• 

<Aq. , 

* 

qj > 

“'^ij/ 1/ j ~3, 

4, •••, (A24) 


We will only consider the eigenvalues and of (Al) 

which coalesce at e = 0. The form of the expansion is 


1/2 


3/2 


a = aj^+e p+ev+e ^ + 


1/2 3/2 

q=q^+E u+ev+e w+ ••* 
If we substitute (A25) in (Al) we obtain 


(A25) 


O(e^); (L - aj_A)qj_ = 0 , 

1/2 

0( e ) : (L - a^A)u = pAqj_ , 


3/2 


(L - Oj^A)v = -Mqj^ + uAu + ( vA + aj^B)qj^ ' 
0( £ ): (L - Oj^A)w = -Mu + uAv + ( vA + aj^B)u 

+ (5A + uB)qj^ . 


(A26 ) 
(A27 ) 
(A28) 


(A29 
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Equation (A26) is the eigenvalue problem for The solution 

of (A27) is u = plus a multiple of the eigenvector qj^ which 

we can take to be zero. The parameter u is still to be 
determined. Next we write 


'' = + nAl + Y. ' 

j'3 


and substitute for v and u in (A28) to obtain 


aifAqi + ^ - tJj_)Aqj = - Mq^_ + u^Aq^^j^ 

j=3 


(A30) 


-f- ( vA + Oj^B)qj^ • 


(A31) 


★ 

If we take the inner product of (A31) with respect to q^^ , 
★ 

and q^ for m = 3, 4, • • • , we obtain, respectively, 


* 


‘^ll 


f 


9 ★ ★ 

p = <Mq^, qj_ > / (a32) 

ail " " ^‘’^'^1' ^11*^ + V + ai<Bqi, qii*> , (^^3) 

^m^'^m - '^l) = ■ q^*> + ai<Bqi, q^*> , 


m = 3, 4, 


(A34) 
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2 

Notice that in order to calculate \i we need only know 
* * 

, and q^^j^ (which is needed in the normalizing condition 
★ 2 

<Aqj^, qj^j^ > = 1). If u is real and y is positive, then the 

perturbed eigenvalues are real Oj^ , .VT Iwl + 0(e), 

2 

and if y is negative then the perturbed eigenvalues are complex 

conjugates 02 = t iVT |y| + O(e^). If y = 0 then also 

u = 0 and the splitting of the eigenvalues occurs at 0(e) or a 

higher order. Assuming that y is real and non-zero, then the 

corresponding eigenvectors of (Al) at this order are q 2 = 

±\fe |yl q^^j^ or q^^, ^2 ^ ^1 “ 1^1 '^ll* Equation (A34) 

provides a linear nonhomogeneous equation for a^j^ and v, while 
the a^ for m = 3, 4, ... are given by (A34). The coefficient a^^ 

ia arbitrary; however if we impose the normalization conditions 
<Aq^^ , q^j^ > = <Aq 2 , qj^j^ > = 1 then a^^ = 0. In order to determine 

V and aj^j^ we turn to (A29) for w. 

We expand w in the series 


w 


biqi . 




(A35) 


substitute in (A29 ) for w, v, and u, and take the inner product 


with respect 

* 

to qj_ . 

We obtain 


II 

0 


qi*> + a^j_ + V + <3qj_^, q^*> 



+ <Bqj_, 

qi*>l- 

(A36) 
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Again, assuming y 0 and solving (A36) and (A33) for v and an we 
obtain 


V 


1 ' 


<Mq 


11 ' 


^1 > 


+ <Mqj,, 


^11 > 





* 


^1 


■I" <Bqj^^ 




(A37) 


an = 1 


<Mq^j_, 







^Bq / 



A <Bq 
2 


1 ' 



(A38) 


This completes the calculation of the eigenvalues and and 

the corresponding eigenvectors through terms 0(e). Note that both 

* ★ 

p and V can be calculated once q^ , qj_ / and q^j^ are known. 

The calculation of higher order corrections for and is 

straightforward. 
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Legends for Figures 

Figure 1. Tne variation of the first five eigenvalues with a 
for 1 < a < 8, n = 0.877 and R = 150. = H 9 . 3 . Eigenvalues 

and are complex conjugates for 2.6 < a < 6.5. The imaginary 

parts of and are shown by the dashed curve with the scale 
given by the axis at the right. The labeling of the eigenvalues 
corresponds to their ordering at R = 0. 

Figure 2. The same as Figure 1 except for n = 0.75 and R = 

85.79 


Figure 3. Tne same as Figure 1 except for n = 0.5 and R = 
68.18 ^ 

Figure 4. The real and imaginary parts of the eigenvalues 0 ^ 
and 0 ^ for n = 0.877, R = 150, and a near a^ = 2.5576. 

O Numerical calculation. Perturbation formula (3.3) 
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Figure 5 Behaviour of the second and third eigenvalues of 
Eqs. (4.2) and (4.3) for n = 0.877 and Ka6"V2 = 0(i). The 
eigenvalues correspond to a double eigenvalue a = -4ir2 for 6 = 0 
and = 0. The perturbation analysis, Eqs. (4.16) and (4.19), is 
given by the dashed curve. The numerical calculation using Eqs. 
(4.2) and (4.3) is shown by the solid curve. On the branch * - ** 
the two eigenvalues are complex conjugates and only the real part 
is shown. 

Figure 6. The behaviour of the first five eigenvalues of Eqs. 
(4.2) and (4.3) for n = 0.877 (5 = 0.123) and 0 < Ra < 1200. The 

eigenvalues and are complex conjugates on the branch * ** 

and only the real part is shown. A branch connects the points O 

and O along which a is real corresponding to two complex 

conjugate values of Ra- 

Figure 7. For points (a, R) between the two curves the 

eigenvalues and of the full problem, Eqs. (1.2)— (1.4), are 

complex conjugates. Numerical calculation. 

Asymptotic formula (4.24). 
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